A Bell pair in a generic random matrix environment 
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Two non-interacting qubits are coupled to an environment. Both coupling and environment are 
represented by random matrix ensembles. The initial state of the pair is a Bell state, though 
we also consider arbitrary pure states. Decoherence of the pair is evaluated analytically in terms 
of purity; Monte Carlo calculations confirm these results and also yield the concurrence of the 
pair. Entanglement within the pair accelerates decoherence. Numerics display the relation between 
concurrence and purity known for Werner states, allowing us to give a formula for concurrence decay. 
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The evolution of entanglement within a pair of qubits 
or spin 1/2 particles under the influence of an environ- 
ment is paradigmatic for the stability of teleportation 
0, and indeed for any quantum information process 0. 
Concurrence provides a measure for the degree of entan- 
glement within such a pair Q . The main purpose of the 
present paper is to establish the generic behavior of the 
decay of entanglement, i.e. concurrence, of such a pair 
of qubits. We propose a random matrix model for the 
unitary evolution of a pair of qubits interacting with one 
or two environments, but not among themselves. The 
environment (s) as well as the couplings will be described 
by one of the classical ensembles of random matrices |^ . 
Research in " quantum chaos" has revealed, that such en- 
sembles describe a chaotic environment well 0, |^ and 
relations to the Caldeira-Legget model have been estab- 
lished 0- In the present article we shall concentrate on 
the Gaussian unitary ensemble (GUE), which describes 
time-reversal breaking dynamics, mainly because it pro- 
vides the simplest analytics. The model we use is based 
on one developed in general for the evolution of deco- 
herence @ and applied to fidelity dec ay El ; the latter 
was successfully tested by experiment [loL llll 0. We 
shall also need information about the evolution of the 
entanglement of the pair with the environment i.e. of 
its decoherence. This we measure in terms of purity [l3| 
rather then von Neumann entropy, because the analytic 
structure of purity allows an analytic treatment in terms 
of a Born expansion. Using this expression we see that 
purity of an entangled state decays faster than purity of 
a product state, but we shall be able to go one step fur- 
ther. Numerically we show that the relation of purity to 
concurrence demonstrated for a specific dynamical model 
[l4| is universal, in the sense that it holds for the random 
matrix model. This relation coincides with the one for a 
Werner state and thus is analytically known, allowing us 
to give a closed, though heuristic, expression for concur- 
rence decay. Both quantities and thus their relation are 
acessable by quantum tomography in experiments with 
trapped ions or atoms, whereinteraction with a controled 



environment is feasable |l5| . 

Concurrence of a density matrix p representing the 
state of a pair of qubits, is defined as 



(1) 



C(p) = max{0, Ai — A2 — A 3 — A 4 } 

where Ai are the eigenvalues of the 
y/ p(o- y <g> Oy)p*(p y (g> Uy) in non- increasing 
(*) denotes complex conjugation in the computational 
basis and a y is a Pauli matrix. Purity is defined as 



matrix 
order; 



P(p) = Tip 2 



(2) 



We study dynamics on a Hilbert space with the struc- 
ture H = Hi ® H\ ® Hi ® H c 2 , where Hf indicates 
(two dimensional) qubit spaces, while H\ will indicate N 
dimensional environments. 

We consider unitary dynamics on the entire space and 
obtain the non-unitary dynamics for the qubits by partial 
tracing over the environment (s). As we wish to consider 
the effect of the environment on the pair of qubits, we 
cannot allow any interaction within the pair, but we con- 
sider interactions with the environments, which may be 
fused to a single one. For convenience we also neglect 
any possible evolution for each qubit individually, which 
is not induced by the coupling to the environment. The 
latter is non-essential to our argument, but simplifies the 
analytic treatment. We thus use the Hamiltonian 



H = H[ + HI + AiU 1 c ' q + X 2 V 2 e ' 



e,q 



(3) 



The first two terms correspond to dynamics of the en- 
vironments. The third and fourth terms represent the 
coupling of each of the qubits to the corresponding en- 
vironment. To obtain further simplification, we consider 
one of the qubits as a spectator, i.e. we assume that it 
has no coupling to an environment (A2 = 0). The corre- 
sponding environment becomes irrelevant and we obtain 
the simplified Hamiltonian 



H x = HI + XVi 



e,q 



(4) 
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Note that we do consider entanglement with the specta- 
tor. This yields the simplest Hamiltonian for which we 
can analyze the effect of an environment on a Bell pair. 
The environment Hamiltonian H° will be chosen from a 
classical ensemble 4] of N x N matrices and the cou- 
pling, Vi' q , from one of 2N x 27V matrices. As usual, the 
GUE ensemble, which represents time-reversal invariancc 
breaking dynamics, is easier to handle analytically than 
the Gaussian orthogonal one. Here we focus on the for- 
mer, while treating the latter in a follow up paper. Evolu- 
tion of both purity and concurrence of the pair of qubits 
can readily be simulated in a Monte Carlo calculation 
and due to the simple structure of purity, it is possible 
to compute analytically this quantity in linear response 
approximation. An exact calculation requires four-point 
functions, which despite of the power of super-symmetric 
techniques [5j are still not readily available. 

To calculate the value of purity, we use the following 
averages and approximations. First we expand the evo- 
lution operator as a Born series; hence we require small A 
and/or short times. We average both V 1 c ' q (which will be 
called V from now on) and HI over the appropriate GUE 
ensemble. Finally we average the initial state and obtain 
Eq. I|12|l . This is the very same scheme followed in || for 
fidelity decay, though details are more complicated Q 
due to the partial traces. 

We define the evolution operator U\ = exp[— iH\t], 
such that the density matrix in Eqs. (Q and is 
p{t) = Tr onv t/ A |V>(0))(V'(0)|t/j, where |^(0)) is the ini- 
tial state of the system. Since Uo is a local operation in 
the environment it will not affect the value of p. Thus 
we can equally evolve with UqU\ instead of U\ alone, ft 
is convenient to use UqU\ since for small A this operator 
will remain in some sense near to unity for longer times. 
We write the Born series to second order as 

Ul U x « 1 - iXI(t) - A 2 J(i); (5) 



I(t) = f AtV(t); J(t) = [ dr / dr'V (t)V(t'). (6) 
Jo Jo Jo 

Here V(t) is the coupling operator in the interaction pic- 
ture: V{t) = U^ VU . Writing |V>(0)) = EJ=i^Im), and 
using Eq. ©, purity reads as: 

P(t)^l-X 2 (RcA J -A 1 -A 2 +ReA 3 ); (7) 

Aj — 4:X fl x*,j k Xi'j'k'X*j, k , Jijk,v(t), (8a) 

Ai = 2xpX„Xi'j/k<x i j, k ,Iijk,ii{t)I i iji s . >I/ (t), (8b) 

A2 = 2x^x*,j k Xi'j'k'X u Iijk i ^(t)I^, k , l ,(t), (8c) 

A 3 = 2x u x*, ]k x v x* J , k ,I ljk ^{t)I l , j , k ,, v (t) (8d) 
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FIG. 1: We show the evolution of purity for separable states 
(purple triangles) , Bell states (green circles) and concurrence 
for Bell states (blue squares) in the crossover regime. The 
lines show theoretical predictions given by Eqs.O and <tHt . 
The environment has dimension 2 10 and the perturbation 
strength is A = 0.025. In the inset we observe the Fermi 
golden rule regime for a larger perturbation A = 0.1. 

(summation over repeated indices is assumed). Indices 
run as follows: Greek ones over the whole Hilbert space, 
the i's over the environment, j's over the first qubit and 
fe's over the spectator qubit. Note that we use the natural 
notation for the indices of vectors in a space which is a 
tensor product of several spaces. 

We now average the perturbation V over the GUE us- 
ing (V m , n ) = and (V m , m 'V n , n >) = 5 m ,n'S m ',n- Due to 
the unitary invariance of the GUE we choose the basis 
that diagonalizes Hf yielding eigenvalues Ei. Then 

\J\tt)) ijk,i' j f k' — *23ijk,i' j' k' 

x f&r f drTe«( T '- T )( E '»- Ei '. (9) 
Jo Jo v , 

The matrix elements of the tensors I ® I and 1 ® I* , 
averaged, yield 

(Iijk,lmn {t\)Ii' j' k ! ,1'm'n' (^2 )) 

= 6--t.„ > 5;-,vi > f dr dr'e i( ' B ^ £;i)(r " r ' ) 

— u ij k,l' m' n u i' 3' k .Iran' j Ll/ / c 
JO JO 

Next we average Hf over the GUE using that 
(Ei,i> = N[l + S(t/r„) - b 2 (t/r H )], where 

&2(i/Tff) is the form factor of the GUE [5j and th the 
Heiscnbcrg time, set to 27T throughout this letter. 

The initial state is a product of pure states for the qubit 
pair and the environment. For the latter we use a random 
initial state X^^il*)) constructed in the large N limit, 
using complex random numbers Xi distributed according 
to a Gaussian centered around zero with width l/\/~N. 
For the pair of qubits we choose a completely general 
pure state. Since we still have the freedom to perform an 
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arbitrary unitary local operation on each qubit, we pick 
a basis in which the initial state for the two qubits is 



\<p a ) = cosa|00) + sina|ll), a <E [0,7r/4]. 



(10) 



The degree of entanglement is characterized by a; in fact 
C{\ip a )(ip a \) — sin2a. Hence our initial state can be 
written as |^(0)) = ^ Xi\i)\(p a ). Neglecting higher order 
terms in 1/N, we obtain ((Aj)) = 2/(t) and ((A 2 )) = 
g a f(t) with 



/(*) = 



2tr h + §£- ifO <t<r H 
2t 2 + ^ if t>T H , 



(11) 



and g a — cos 4 a + sin 4 a. To leading order ((Ai)) = 
((A 3 )} = 0. We obtain 



PMt) = \-\ 2 {2-g a )f{t). 



(12) 



From this result we see directly that purity decay will 
be faster the more entangled the initial state was. The 
validity of this approximation is limited to large values 
of purity, i.e. short times or weak coupling. This is 
valuable for applications to quantum information, but 
we are interested in the dynamical picture as a whole 
and thus would like to obtain an expression valid for a 
wide range of physical situations. As a way to achieve 
this for fidelity decay, exponentiation of the leading term 
of the linear response formula was proposed |l(j | . A sim- 
ilar approximation is taken here, to obtain Pelr^)- In 
order to calculate the appropriate formula, we must sat- 
isfy Pelr(£) ~ Pis.it) for small t, and consider correct 
asymptotics. These will be estimated as the purity af- 
ter applying a totally depolarizing channel on one qubit 
to the 2 qubit state Eq. lJTf)|l. The expected asymptotic 
value is g a /2, and the final expression is 



lxr(*) = 



i _ i^. 

2 



e i-so/2 



(13) 



This result is in excellent agreement with numerics as 
shown in Fig. ^ and displays the transition from expo- 
nential to Gaussian decay as Heisenberg time is reached. 

We thus have an approximate formula for the decay 
of purity of a Bell pair. What about concurrence? At 
this point we take up a result [l4| for the behavior of a 
Bell pair coupled to a kicked spin chain 01 • F° r a wide 
range of situations the decay of a pure Bell state leads to 
purities and concurrences that closely follow those of a 
Werner state in a Concurrence-Purity (CP) diagram [lif . 
To test model independence, and thus universality of this 
behavior we make the corresponding numerical simula- 
tions in the RMT model. We find that the Werner state 
CP relation is quite well fulfilled in the large N limit, 
as can be seen in Fig. [21 where results for fixed coupling 
but different sizes of the RMT environment are shown. 
Studying other couplings allowed by the full Hamiltonian 
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FIG. 2: The relation between C and P for a fixed coupling 
(A = 0.3) and variable size of the environment. We average 
over 10 realizations of the Hamiltonian Q and 15 initial con- 
ditions. The gray area indicates the region of physical states, 
and the line the concurrence-purity relation for Werner states. 
In the inset we observe short time deviations from this rela- 
tion. 



Eq. (|3J) leads to similar results even if we are in differ- 
ent purity decay regimes. A partial explanation for this 
behavior can be found in |l8| . 

We thus have the second relevant result of this pa- 
per; namely the relation of purity to concurrence for a 
non-interacting Bell pair interacting with a chaotic envi- 
ronment follows generically the curve of a Werner state. 
The importance of this statement is underlined by the 
fact, that the actual state reached at any time is typi- 
cally not a Werner state. This is tested, by considering 
the spectrum of the density matrix, which should display 
a triple degeneracy for a Werner state. In fact a typical 
spectrum at P = 0.51 is {0.692, 0.142, 0.110, 0.056}. 

Having established the genericity of the Werner state 
relation we can now insert the expression (|13fl for purity 
into the latter and obtain the heuristic expression 



C E lr(<) = max < 0, 



y/l2P ELR (t) - 3 - 1 



(14) 



for concurrence decay. In Fig. [3] we can see that this re- 
lation is well obeyed by Monte Carlo calculations. Here 
we obtain two different time regimes, an exponential one 
(Fermi golden rule) for strong perturbations and a Gaus- 
sian one for weak perturbation. The time scale which de- 
fines the crossover between the two regimes is the Heisen- 
berg time of the environment. Since the exponential be- 
havior can be obtained letting the Heisenberg time go to 
infinity, we retrieve results derived from a master equa- 
tion approach (lij . 

We have thus obtained a satisfactory expression for 
concurrence decay, but we have to remember that we 
used an extremely simplified model. In some points the 
linear response treatment of the RMT model is slightly 
affected by these assumptions; we have made calculations 



4 




5 10 15 20 



t/T H 

FIG. 3: The average evolution of concurrence using and en- 
vironment of dimension 128 and averaging over 10 Hamilto- 
nians and 15 initial conditions is shown. The black curve 
corresponds to Eq. 11411 . We show both the Gaussian regime 
(S = 0.008) and the Fermi golden rule regime in the inset 
(6 = 0.07). 

where these approximations where lifted. In other words 
we have allowed both qubits to interact, either with inde- 
pendent environments or with the same one, see Eq. 
The form of our result is essentially the same: 

PLr® = 1 - (2 - g a )[\lh(t) + Al/ 2 (t)], (15) 

where fi(t) is identical to f(t) as in Eq. but us- 

ing the Heisenberg time of Hf. Other generalizations are 
possible. Local Hamiltonians for each qubit causing the 
degeneracy of the levels of the qubits to be lifted can be 
included. Environment Hamiltonians and couplings cho- 
sen from the Gaussian orthogonal ensemble may also be 
considered, as well as mixed states for the environment. 
Detailed linear response calculations as well as Monte 
Carlo calculations for these cases will be presented in the 
follow up paper. 

Some points are worth mentioning: a) We see signifi- 
cant deviations from the usual exponential decay at times 
of the order of the Heisenberg time as defined by the en- 
vironment. Thus, if the spectrum of the environment 
becomes very dense, and correspondingly the Heisenberg 
time moves off to infinity we recover the usual stochastic 
result, b) If the transition region can actually be seen, 
then the spectral stiffness of a chaotic environment has 
a small but significant stabilizing effect. The absence of 
spectral stiffness can be modeled by the so-called Poisson 
random ensemble [2(| • c) We have limited our discussions 
to the GUE for two reasons. The simple form of the form 
factor yields a concise final expression for purity decay. 
An additional advantage resulting from the unitary in- 
variance of the coupling term, is that the final result is 
invariant under any local operation at each qubit. This 
is no longer guaranteed for orthogonal invariance only. 



The implications will be discussed in another paper. 

Summarizing, we have developed a random matrix 
model for the evolution of a Bell pair interacting with 
a generic chaotic environment. Within this model we 
derive the linear response approximation for the purity 
decay of a Bell pair and show that it differs significantly 
from decay of a product state of two spins or qubits, 
even in the extreme case, where one of the qubits is only 
a spectator. Exponentiation extends the validity of this 
result far beyond its original reach. Monte Carlo calcu- 
lations show that the relation between concurrence and 
purity, as obtained for Werner states, holds for RMT 
models and we thus expect it to be generic. Based on 
these results we have obtained and tested a heuristic for- 
mula for the decay of concurrence of a non-interacting 
Bell pair. 
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